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I. INTRODUCTION 

Dissociation, fragmentation, and fissioning processes underly physical and chemical phe- 
nomena in a variety of finite-size systems, characterized by a wide spectrum of energy scales, 
nature of interactions, and characteristic spatial and temporal scales. These include nu- 
clear fission unimolecular decay and reactions in atoms and molecules 0, and more 
recently dissociation and fragmentation processes in atomic and molecular clusters 
Investigations of the energetics, mechanisms, pathways, and dynamics of fragmentation 
processes provide ways and means for explorations of the structure, stability, excitations, 
and dynamics in the many-body finite systems mentioned above, as well as they allow for 
comprehensive tests of theoretical methodologies and conceptual developments, and have 
formed active areas of fruitful research endeavors in nuclear physics, and more recently in 
cluster science. 

Under the general title of dissociation and fragmentation [0 processes in metal clus- 
ters, one usually distinguishes two classes of phenomena, i.e., (1) dissociation of neutral 
monomers and/or dimers, and (2) fission. The physical processes in the first class are most 
often referrred to as evaporation of monomers and/or dimers, since they are endothermic 
processes and are usually induced through laser heating of the cluster. The unimolecular 
equations associated with these processes are 

M+ M+_i + M , (1) 



for monomer separation, and 



M+ ^ M+_2 + M2 , (2) 



for dimer separation (A^ denotes the number of atoms in the clusters 0). The parent 
clusters have been taken here as being singly ionized, in order to conform with available 
experimental measurements |Q. Fission on the other hand, is most often an exothermic 
process and is due to the Coulombic forces associated with excess charges on the cluster. It 



* Contributed Chapter to the book Metal Clusters, Edited by W. Ekardt (Wiley, New York, 
1999) pp. 145-180. 
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has been found that the minimum excess charge required to induce fission is 2 elementary 
units (either positive or negative). In this case the doubly-charged parent cluster splits into 
two singly charged fragments, and the corresponding unimolecular equation can be written 
as 

+ p, P = 1, [N/2] . (3) 

It needs to be emphasized that fragmentation through fission involves most often the over- 
coming of a fission barrier, while momomer and dimer separation are barrierless processes 

i- 



A. Metal cluster fission and nuclear fission: Similarities and differences 

Multiply charged metallic clusters (M^"*") are observable in mass spectra if they exceed 
a critical size of stability N^^ (e.g. for Z = 2, Nl+ = 27 for Na and = 20 for K 
For clusters with > N^~^, evaporation of neutral species is the preferred dissociation 
channel, while, below the critical size, fission into two charged fragments dominates (for 
Z = 2, two singly charged fragments emerge). Nevertheless, at low enough temperature, 
such M|+ {N < iVf +) clust ers can be metastable above a certain size N^~^, because of the 
existence of a fission barrier Ef, (for Na^ and K^, A*"^"*' = 7 [p!0| , pJ]] ) . 

These observations indicate that fission of metal clusters occurs when the repulsive 
Coulomb forces due to the accumulation of the excess charges overcome the electronic 
binding (cohesion) of the cluster. This reminds us immediately of the well-studied nuclear 
fission phenomenon and the celebrated Liquid Drop Model (LDM) according to which the 
binding nuclear forces are expressed as a sum of volume and surface terms, and the balance 
between the Coulomb repulsion and the increase in surface area upon volume conserving 
deformations allows for an estimate of the stability and fissility of the nucleus [ P^ , P^ . 

We note that for doubly charged metal clusters with < 12 microscopic descriptions of 
energetics and dynamics of fission, based on first-principles electronic-structure calculations 
in conjunction with molecular dynamics (MD) simulations, have been performed [^jOI (^^^ 
section III.C.l for details). Several of the trends exhibited by the microscopic calculations 
(such as infiuence of magic numbers, associated with electronic shell closing, on fission 
energetics and barrier heights; predominance of an asymmetric fission channel; double- 
humped fission-barrier shapes; shapes of deforming clusters along the fission trajectory 
portraying two fragments connected through a stretching neck) suggest that appropriate 
adaptation of methodologies developed originally in the context of nuclear fission may 
provide a useful conceptual and calculational framework for studies of systematics and 
patterns of fission processes in metallic clusters. 

In this context, it is useful to comment on the earliest treatments of pertinent nuclear 
processes, i.e., fission and alpha radioactivity |jl^,|l5|,||. Adaptation of the simple one- 



center LDM to charged metallic clusters , involving calculation of the Coulomb repulsive 
energy due to an excess charge localized at the surface, yields a reduced LDM fissility 
parameter ^ = {Z'^/N)/{Z'^/N)cr, where {Z'^/N)cr = IGTrrfcr/e^ with the surface energy 
per unit area denoted by a and being the Wigner-Seitz radius (using bulk and a 
values, {Z'^/N)cr = 0.44 and 0.39 for K^"*" and Na^"*", respectively). Accordingly, a cluster 
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is unstable for > 1 (implying that for with < 9 and Na^ with < 10 barrierless 
fission should occur) with the most favorable channel being the symmetric one (i.e., when 
the two fragments have equal masses, which is only approximately true for nuclear fission, 
and certainly not the case for small metal clusters). For 0.351 < ^ < 1, the system is 
metastable (i.e., may fission in a process involving a barrier), and for < ^ < 0.351 the 
system is stable. 

At the other limit, a-radioactivity, which may be viewed as an extreme case of (su- 
perasymmetric) fission, is commonly described as a process where the fragments are formed 
(or as often said, preformed) before the system reaches the top of the barrier (saddle point), 
and as a result the barrier is mainly Coulombic . We note here that asymmetric emission 
of heavier nuclei is also known (e.g., ^^^Ra— >^^C+^°^Pb, referred to as exotic or cluster 
radioactivity |[T6|-[18|), and the barriers in these cases resemble the one- humped barrier of 
alpha radioactivity and do not exhibit modulations due to shell effects |l^ . We also remark 
that such a-radioactivity-type (essentially Coulombic) barriers have been proposed recently 
19| for describing the overall shape of the fission barriers in the case of metal clusters. 

Although, several aspects of the simple LDM (e.g., competition between Coulomb and 
surface terms) and the a-particle, Coulombic model (e.g., asymmetric channels and a scis- 
sion configuration close to the location of the saddle of the multi-dimensional potential- 
energy surface) are present in the fission of metal clusters, neither model is adequate in light 
of the characteristic behavior revealed from the microscopic calculations and experiments. 
Rather, we find that proper treatments of fission in these systems require consideration of 
shell effects (for a recent experimental study that demonstrates the importance of shell ef- 
fects in metal-cluster fission, see Ref. [9b]). While such effects are known to have important 
consequences in nuclear fission (transforming the one-humped LDM barrier for symmetric 
fission into a two- humped barrier pO|J^), their role in the case of metal clusters goes even 
further. Indeed, as illustrated below (see section II1.C.2) for the case of the magic Na^g (8 
delocalized electrons), shell effects can be the largest contribution to the fission barrier, in 
particular in instances when the LDM component exhibits no barrier (in this case the LDM 
fissility ^ > 1). In this respect, Na^Q is analogous to the case of superheavy nuclei, which 
are believed pl| to be stabilized by the shell structure of a major shell closure at Zp = 114, 
Nn = 184 {Zp is the number of protons and Nn is the number of neutrons; unfortunately 
such nuclei have not been yet observed or synthesized artificially). 

B. Other decay modes in atomic and molecular clusters 

In this chapter, we will concentrate on the unimolecular processes in metal clusters 
described by Eqs. (|l|— |). However, there is a variety of additional dissociation and frag- 



mentation modes in atomic and molecular clusters (see reviews in Ref. [22 ), which have 



been discovered experimentally or anticipated theoretically; among them we mention: 
1. Unimolecular fission of triply and higher charged cationic simple metal clusters 



2. Metastability against electron autodetachment of multiply charged anionic atomic 
clusters [p5|-p^ and fuUerenes [p6|- p^] ; 
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3. Fragmentation of cationic fullerenes via sequential evaporation of carbon dimers p9 



4. Ultrarapid fragmentation of rare-gas clusters following excitation (involving excimer 



formation 30 ) or ionization 22 



5. Multifragmentation phase transitions according to microcanonical thermodynamics 
of highly excited atomic clusters ; and 

6. Pathways and dynamics of dissociation and fragmentation of ionized Van-der-Waals 



and hydrogen-bonded molecular clusters p2|,|32 



C. Organization of the chapter 



In the following, we will present jellium-related theoretical approaches [specifically the 
Shell Correction Method (SCM) and variants thereof] appropriate for describing shell ef- 
fects, energetics and decay pathways of metal-cluster fragmentation processes (both the 
monomer/dimer dissociation and fission), which were inspired by the many similarities 
with the physics of shell effects in atomic nuclei (section II). In section III, we will compare 
the experimental trends with the resulting theoretical SCM interpretations, and in addi- 
tion we will discuss theoretical results from first-principles MD simulations (section III.C.l). 
Section IV will discuss some latest insights concerning the importance of electronic-entropy 
and finite-temperature effects. Finally, section V will provide a summary. 



II. THEORY OF SHAPE DEFORMATIONS 

In early applications of the jellium model, the shape of metal clusters was assumed in all 
instances to be spherical P5|J5^ , but soon it became apparent that the spherical symmetry 
was too restrictive [^,^]. Indeed clusters with open electronic shells (between the magic 



numbers N^. = 2, 8, 20, 40, 58, 92, etc..) are subjected to Jahn- Teller distortions |^ . 
By now it has been well established that a quantitative description of the underlying shell 
effects and of fragmentation phenomena (as well as of other less complicated phenomena 
such as Ionization and Vertical Electron Detachment) requires a proper description of the 
deformed shapes of both parent and daughter clusters (of both precursor and final ionic or 
neutral product in the case of ionization and vertical electron detachment). 

A most successful method for describing both deformation and shell effects in simple 
metal clusters (i.e., those that can be described by the jellium background model) is the 
SCM, originally developed in the field of nuclear physics |28|,||. In a series of recent publica- 
tions p5| , p6| , p8| , |39| -^ , the SCM was further developed, adapted, and applied in the realm 
of finite-size, condensed- matter nanostructures (i.e., metal clusters [p5|, p6| , p^ -^ , but also 
multiply charged fullerenes [^, ^He clusters |^4[, and metallic nanowires and nanoconstric- 



tions |J3|). Additionally, Refs. [^^-^ have used semiempirical versions (see below) of the 
SCM to study the shapes of neutral Na clusters ||46|j47[] and aspects of metal-cluster fission 

El®- 



The SCM derives its justification from the local-density-approximation (LDA) functional 
theory and has been developed as a two-level method. 
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At the microscopic level, referred to as the LDA-SCM, the method has been shown to 
be a non-selfconsistent approximation to the Kohn-Sham (KS) — LDA approach [^. Apart 
from computational efficiency, an important physical insight provided by the LDA-SCM is 
that the total KS-LDA energy Etota\{N) [or in another notation i^^Ks(^)] of a finite system 
of interacting delocalized electrons (or more generally of other fermions, like nucleons or 
^He atoms) can be divided into two contributions, i.e., 

Etotai(A^) = ^(A^) + AE,^{N) , (4) 

where E is the part that varies smoothly as a function of the system size (e.g., the num- 
ber, N of atoms in a metal cluster), while AEsh{N) is an oscillatory term accounting for 
the shell effects; it arises from the discretization of the electronic states (quantum size ef- 
fect). AEsh{N) is usually called a shell correction in the nuclear P^JI[| and cluster pB| , pB 
literature. 

Starting from the fundamental microscopic separation in Eq. (^), various semiempirical 
implementations (referred to as SE-SCM, see section II. B) of such a division consist of 
different approximate choices and methods for evaluating the two terms contributing to 
this separation. 

As an illustration of the physical content of Eq. (|^) (which as well serves as a motivating 
example for the SCM), we show in Fig. 1 the size-evolutionary pattern of the Ionization Po- 
tentials (IPs) of NaAT clusters, which exhibits odd-even oscillations in the observed spectrum 
in addition to the major features (major IP drops) at the magic numbers. Theoretical cal- 
culations at three different levels are contrasted to the experimental observations, namely, 
a smooth description of the pattern [Inset (a)], and two levels of shell-corrected descriptions 
— one assuming spherical symmetry [Inset (b)], and the other allowing for triaxial shape 
deformations [Fig. 1, main frame]. The progressive improvement of the level of agreement 
between the experimental |51|J5^ and theoretical patterns is evident. 

Below, we first outline the microscopic derivation of Eq. (^, and subsequently we pro- 
ceed with a presentation of the SE-SCM. 



A. Microscopic Foundation of Shell Correction Methods The LDA-SCM 

The LDA-SCM approach, which has been shown to yield results in excellent agreement 
with self-consistent KS-LDA calculations [^,^], is equivalent to a Harris functional ||53 



approximation (-EHarris[p"], see below) to the KS-LDA total energy [^ (EksIprs]), with the 
input density obtained through a variational minimization of an extended Thomas- Fermi 
(ETF) energy functional, Eetf[p]- 

The property of the non-selfconsistent Harris functional to yield total energies close to 
the KS-LDA ones is based on the following equality: 

^KS[PKS] = Euarrisir] + 0{6p^) , (5) 

where Sp = pks — P™- Namely, the KS-LDA energy is, to second-order in 5p, equal to the 
Harris energy. 



Several recent publications have proven [54-56] the validity of equation (B[) in connec- 



tion with the Harris functional, which is often used in electronic structure calculations of 
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FIG. 1. IPs of NaAT clusters. Open squares: Experimental measurements [^,52]. Solid 
circles: Theoretical IPs derived from the SCM assuming ellipsoidal (triaxial) deformations. Inset 
(a): The solid line represents the smooth contribution to the theoretical SCM IPs. Inset (b): The 
solid circles are the IPs derived from the SCM assuming spherical symmetry. 

molecules, surfaces, and other condensed- matter systems. We note that, in the context 
of nuclear physics, Strutinsky had earlier proven ||3^ the validity of Eq. (|^), with the dif- 
ference that he utilized the Hartree-Fock (HF) functional instead of the KS-LDA one. In 
the nuclear-physics literature, the HF version of Eq. (^) is referred to as the Strutinsky 
theorem. 

Usually, in the Harris functional, the input density is taken as a superposition of 
site densities. Initially the site components of the input density were not optimized. 
Later [0,^], it was realized that the results could be improved by variationally adjusting 



the site components through a maximization of the Harris functional itself. However, doing 
so adds the burden of a matrix diagonalization for obtaining the eigenvalues (see below) 
at each step of the variation. Our method differs from the Harris approach in that the 
optimization of the input density is achieved by us through a variational ETF method |57 
(which does not require such a step- by-step matrix diagonalization). 

The non-selfconsistent Harris functional is given by the following expression. 



-E' Harris [P™] — 

+ E - / {Ivh[p"\t)] + Kc[p^°(r)]} p''^(r)t/r + J S^,[p"\r)]dr , (6) 



where Vh is the Hartree (electronic) repulsive potential, Ei is the repulsive electrostatic 
energy of the ions, and E^dp] = J ^xc[p]dr is the exchange-correlation (xc) functional pS 



[the corresponding xc potential is given as V^c(r) = ^-E'xc[p]/'5p(r)]. e°"* are the eigenvalues 
(non-selfconsistent) of the single-particle Hamiltonian, 



H 



2m, 



with the mean-field potential given by 



(7) 



(8) 



Vi{r) being the attractive potential between the electrons and ions. 

The ETF-LDA energy functional, E-etf[p], is obtained by replacing the kinetic energy 
term, T[p\, in the usual LDA functional, namely in the expression. 



-EldaH = 

T[p] + / [^VH[pir)] + Vjir)^pir)dv + J £^Mr)]dr + Ej 



(9) 



by the ETF kinetic energy, given to the 4th-order gradients as follows [pQ 

tETF[p]dr 

1 (Vp)2 



TetfIp] 



2mp 



1 Yl 

3 V P . 



(3^2)2/3^5/3 + 

-3d p 



1 

270 

2 



(Stt 
dv 



2a-2/3^1/3 



(10) 



We would like to remind the reader that the KS kinetic energy is of course given by the 
expression 



Tks[p 



KS 



E< 



?KS,i 



2mp 
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where the single-particle wave functions 0KS,i(r) are obtained from a self-consistent solution 
of the KS equations. 

The optimal ETF-LDA total energy is obtained by minimization of -Eetf [p] with respect 
to the density. In our calculations, we use for the trial densities parametrized profiles 
p(r; {7i}) [l60| , p5i^ with {■ji} as variational parameters (the ETF-LDA optimal density 
is denoted as p). The single-particle eigenvalues, {e°"*}, in Eq. (||) are obtained then as 
the solutions to the single-particle Hamiltonian of Eq. (|^) with V^n replaced by Vetf [given 
by Eq. (^ with p™(r) replaced by p(r)]. Hereafter, these single-particle eigenvalues will be 
denoted by {et}. 

In our approach, the smooth contribution in the separation (^) of the total energy is 
given by Eetf[p], while the shell correction, AEsh_, is simply the difference p5|j26[1 

A-E'sh= Eua.ms[p] — -E'eTf[p] 
occ „ 

= E - y p(r)V^ETF(r)rfr - Tetf[p] • (12) 



i=l 
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B. Semiempirical shell-correction method (SE-SCM) 



1. Methodology 

Rather than proceed with the microscopic route, Strutinsky proposed a method for the 
separation of the total energy into smooth and shell- correct ion terms [see Eq. (^)] based on 
an averaging procedure. Accordingly, a smooth part, Egp, is extracted out of the sum of 
the single-particle energies J^T^^i [or S^'^^i"*' see Eq. (^] by averaging them through an 
appropriate procedure. Usually, but not necessarily, one replaces the delta functions in the 
single-particle density of states by gaussians or other appropriate weighting functions. As 
a result, each single-particle level is assigned an averaging occupation number /j, and the 
smooth part Esp is formally written as 

Esp = J2eJi. (13) 

i 

Consequently, the Strutinsky shell correction is given by 

occ 

Ai?i*^ = E^^-^sp. (14) 
1=1 



The Strutinsky prescription (|lj) has the practical advantage of using only the single- 
particle energies ej, and not the smooth density p. Taking advantage of this, the single- 
particle energies can be taken as those of an external potential that empirically approx- 
imates the self-consistent potential of a finite system. In the nuclear case, a modified 
anisotropic three-dimensional harmonic oscillator has been used successfully to describe 
the shell-corrections in deformed nuclei [0,0. 

The single-particle smooth part, Egp, however, is only one component in the smooth 
contribution E[p], which needs to be added to the shell correction term in order to yield 
the total energy, i.e., 

Etotal- AES*=^ + E[p] . (15) 

Strutinsky did not address the question of how to calculate microscopically the smooth 
part E (which necessarily entails specifying the smooth density p). Instead he circumvented 
this question by substituting for E the empirical energies, -Eldm, of the nuclear liquid drop 
model, namely he suggested that 

Etotal ~ AE^^' + Eldm . (16) 

In applications of Eq. (0), the single-particle energies involved in the averaging [see 
Eqs. (|l^) and (0)] are commonly obtained as solutions of a Schrodinger equation with 
phenomenological one-body potentials. This last approximation has been very successful 
in describing fission barriers and properties of strongly deformed nuclei using harmonic- 
oscillator-type or Wood-Saxon empirical potentials. 
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2. Liquid- drop model for neutral and charged deformed clusters 



For neutral clusters, the LDM expresses [^,^,^ the smooth part, E, of the total energy 



as the sum of three contributions, namely a volume, a surface, and a curvature term, i.e., 

E £'vol ~l~ -^surf ~l~ -^curv 

= Ay J dr + a J dS + Ac J dSK , (17) 

where dr is the volume element and dS is the surface differential element. The local 
curvature k is defined by the expression k = 0.5 (i?"^^ + R^l^), where -Rmax and -Rmin are 
the two principal radii of curvature at a local point on the surface of the jellium droplet (of a 
general shape) which models the cluster. The corresponding coefficients can be determined 
2^,|2B|BO] by fitting the ETF-LDA total energy for spherical shapes (see section II. A) to the 



following parametrized expression as a function of the number, N, of atoms in the cluster 
0, 



The following expressions relate |^ the coefficients cr, and A^ to the corresponding 
coefficients, (a's), in Eq. (p!^), 

3 11 

Av = - — ra^, ; cr = - — -as ; Ac = Oc ■ (19) 

47rrf Airrf AirVs 

In the following, we will focus on the case of clusters with ellipsoidal (triaxial) shapes. 
In the case of ellipsoidal shapes the areal integral and the integrated curvature can be 
expressed in closed analytical form with the help of the incomplete elliptic integrals J-'i'ip, k) 
and S{ijj, k) of the first and second kind respectively. Before writing the formulas, we 
need to introduce some notations. Volume conservation must be employed, namely 

a'h'c'/Rl = a6c = 1 , (20) 

where Rq is the radius of a sphere with the same volume {Rq = VsN^^'^ is taken to be 
the radius of the positive jellium assuming spherical symmetry, being the corresponding 
Wigner-Seitz radius), and a = a'/RQ, etc., are the dimensionless semi-axes. The eccen- 
tricities are defined through the dimensionless semi-axes as follows 

el = l- {c/af 
el = l- {h/af 

el = l- {c/hf . (21) 



The semi-axes are chosen so that 



a>b>c. (22) 



With the notation sinip = ei, k2 = 62/61, and k^ = 63/61, the relative (with respect to 
the spherical shape) surface and curvature energies are given by 
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and 



-^surf 
T^sph 
-^surf 



<rv _ be 
Ecurv 2(3 



ab 
~2 



l-e| 

ei 



J^{ij,h) + ei£{ij,h)+e' 



(23) 



1 + - ((1 - e?)^(^, fcs) + e'iS{^P, h 
^1 



(24) 



The change in the smooth part of the cluster total energy due to the excess charge ±Z 
has been discussed for spherical clusters in Refs. p5| , p6i . The result may be summarized as 



AE^p^(z) = s^p'^(z) - £;^p'xo) = twz + 



Z(Z±0.25)e2 
2(i?o + 5) 



(25) 



where the upper and lower signs correspond to negatively and positively charged states, 
respectively, W is the work function of the metal, Rq is the radius of the positive jellium 
assuming spherical symmetry, and 5 is a spillout-type parameter. 

To generalize the above results to an ellipsoidal shape, 0(-Ro + 5) = e^ /{Rq + 5), which 
is the value of the potential on the surface of a spherical conductor, needs to be replaced 
by the corresponding expression for the potential on the surface of a conducting ellipsoid. 
The final result, normalized to the spherical shape, is given by the expression 



AE^"(Z) ± WZ 
AE^p^{Z) ± WZ 



ei 



(26) 



where the ± sign in front of WZ corresponds to negatively and positively charged clusters, 
respectively. 



3. The modified Nilsson potential for ellipsoidal shapes 

A natural choice for an external potential to be used for calculating shell corrections 
with the Strutinsky method is an anisotropic, three-dimensional oscillator with an addi- 
tional F angular-momentum term for lifting the harmonic oscillator degeneracies |6^. Such 
an oscillator model for approximating the total energies of metal clusters, but without sepa- 
rating them into a smooth and a shell-correction part in the spirit of Strutinsky' s approach. 



had been used |52] with some success for calculating relative energy surfaces and deforma- 



tion shapes of metal clusters. However, this simple harmonic oscillator model had serious 
limitations, since i) the total energies were calculated by the expression fZ^^^'^ei, and thus 
did not compare with the total energies obtained from the KS-LDA approach, and ii) the 
model could not be extended to the case of charged (cationic or anionic) clusters. Thus 
absolute ionization potentials, electron affinities, and fission energetics could not be calcu- 
lated in this model. Alternatively, in our approach, we are making only a limited use of 
the external oscillator potential in calculating a modified Strutinsky shell correction. Total 
energies are evaluated by adding this shell correction to the smooth LDM energies (which 
incorporate xc contributions, since the LDM coefficients are extracted via a comparison 
with total ETF-LDA energies, or they are taken from experimental values). 
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In particular, a modified Nilsson Hamiltonian appropriate for metal clusters [p5| , p6 
given by 



IS 



(27) 



where Hq is the hamiltonian for a three-dimensional anisotropic oscillator, namely 



Y^{alak + ]-)huj,. (28) 

fc=i 



Uq in Eq. (|27|) is a dimensionless parameter, which for occupied states may depend on 
the effective principal quantum number n = ni + n2 + ^3 associated with the major shells 
of any spherical-oscillator, (?T,i,n2,n3) being the quantum numbers specifying the single- 
particle levels of the hamiltonian Hq (for clusters comprising up to 100 valence electrons, 
only a weak dependence on n is found, see Table I in Ref. [40a]). Uq vanishes for values of 
n higher than the corresponding value of the last partially (or fully) filled major shell with 
reference to the spherical limit. 

]? = Y,k=i is ^ " stretched" angular momentum which scales to the ellipsoidal shape 
and is defined as follows, 

ll = iqiP2-q2Pi)\ (29) 



(with similarly obtained expressions for li and I2 via a cyclic permutation of indices) where 
the stretched position and momentum coordinates are defined via the corresponding natural 
coordinates, g^*^* and p^^*, as follows, 

qk ^ qT\rneUJk/hfl^ = , (A: = 1, 2, 3) , (30) 



p, ^ pni/hm^u,)'/' = , (fc = 1, 2, 3) . (31) 

The stretched F is not a properly defined angular-momentum operator, but has the 
advantageous property that it does not mix deformed states which correspond to sherical 
major shells with different principal quantum numbers n = ni + n2 + n3 (see, the Appendix 
in Ref. [40a] for the expression of the matrix elements of F). 

The subtraction of the term < P >«= n(n^ where < denotes the expectation 
value taken over the nt/i-major shell in spherical symmetry, guaranties that the average 
separation between major oscillator shells is not affected as a result of the lifting of the 
degeneracy. 

The oscillator frequencies can be related to the principal semi-axes a', 6', and d [see, 
Eq. (^)] via the volume-conservation constraint and the requirement that the surface of 
the cluster is an equipotential one, namely 
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ujia' = uj2b' = ujsc = ujqRo , 



(32) 



where the frequency ujq for the spherical shape (with radius Rq) was taken according to 



Ref. to be 

49 eV bohr^ 



r.2jYl/3 



1+ ' 



r 



Ari/3 



(33) 



Since in this paper we consider solely monovalent elements, in Eq. (|3^) is the number 
of atoms for the family of clusters M^"^, is the Wigner-Seitz radius expressed in atomic 



units, and t denotes the electronic spillout for the neutral cluster according to Ref. |]35[| . 



4- Shell correction and averaging of single-particle spectra for the modified Nilsson potential 



Usually Egp [see Eqs. (|T3D and (|T^] is calculated numerically |Q. However, a variation 
of the numerical Strutinsky averaging method consists in using the semiclassical parti- 
tion function and in expanding it in powers of h . With this method, for the case of an 
anisotropic, fully triaxial oscillator, one finds [0,^ an analytical result, namely ||68| 



if3Ar ^4/3 , 1 ujI + + ujI 2/3' 
where Np denotes the number of delocalized valence electrons in the cluster. 



- (i(3iV.)'" + ^^L±^(3iV.)=/3 ) . (34) 



In the present work, expression (|34D (as modified below) will be substituted for the 
average part Esp in Eq- (P^, while the sum J^T^^i will be calculated numerically by spec- 
ifying the occupied single-particle states of the modified Nilsson oscillator represented by 
the hamiltonian (P?]). 

In the case of an isotropic oscillator, not only the smooth contribution, E°^'^, but also 
the Strutinsky shell correction (|T^) can be specified analytically pi] with the result 



A^sto(^) = ^^^o(3iV,)2/=^(-l + 12x(l - x)) , (35) 

where x is the fractional filling of the highest partially filled harmonic-oscillator major 
shell. For a filled shell (x = or 1), AEf!o(0) = ~M^^oi^N^)^^^, instead of the essentially 
vanishing value as in the case of the ETF-LDA defined shell correction (cf. Fig. 1 of Ref. 
[40a]). To adjust for this discrepancy, we add — Aii^fiJ^Q(O) to AEf-^^ calculated through Eq. 
(|l^) for the case of open-shell, as well as closed-shell clusters. 



5. Overall procedure 

We are now in a position to summarize the calculational procedure for the SE-SCM in 
the case of deformed clusters, which consists of the following steps: 

1. Parametrize results of ETF-LDA calculations for spherical neutral jellia according to 
Eq. (pp. 
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2. Use above parametrization (assuming that parameters per differential element of vol- 
ume, surface, and integrated curvature are shape independent) in Eq. (|l^) to calculate 
the hquid-drop energy associated with neutral clusters, and then add to it the charg- 
ing energy according to Eq. ( ^61) to determine the total LDM energy E (available 
experimental values for a and W can also be used). 

3. Use Equations (^) and (^Sf ) for a given deformation [i.e., a', b', c', or equivalently ui, 
LJ2, t^3, see Eq. (|3^)] to solve for the single-particle spectrum (e^). 

4. Evaluate the average, Egp, of the single-particle spectrum according to Eq. ( p4D and 
subsequent remarks. 

5. Use the results of steps 3 and 4 above to calculate the shell correction AEf^^ according 
to Eq. dll). 

6. Finally, calculate the total energy -E'totai as the sum of the liquid-drop contribution 
(step 2) and the shell correction (step 5), namely -Etotai = E + AE^-^^. 

The optimal ellipsoidal geometries for a given cluster M^^, either neutral or charged, 
are determined by systematically varying the distortion (namely, the parameters a and b) in 
order to locate the global minimum of the total energy ii^totai(^; Z) (for the global minima 
and equilibrium shapes of neutral NaAr clusters according to the ellipsoidal model in the 
range 3 < < 60, see Fig. 22 of Ref. [40a]). 



6. Asymmetric two-center oscillator model for fission 

Naturally, the one-center Nilsson oscillator is not the most appropriate empirical po- 
tential for describing binary fission, which involves the gradual emergence of two separate 
fragments. A better choice is the asymmetric two-center oscillator model (ATCOM). Ac- 
cording to the ATCOM approach, the single-particle levels, associated with both the initial 
one-fragment parent and the separated daughters emerging from binary cluster fission, are 



determined by the following single-particle hamiltonian ||69| , |70 



H = T + ^meuj%p^ + \m,ujllz - Zif + Kcck(2) + f/(l?) , (36) 

where i = \ for 2; < (left) and i = 2 for 2; > (right). 

This hamiltonian is axially symmetric along the z axis, p denotes the cylindrical coor- 
dinate perpendicular to the symmetry axis |7l[|. The shapes described by this Hamiltonian 
are those of two semispheroids (either prolate or oblate) connected by a smooth neck [which 
is specified by the term Kieck(^)]- -Zi < and ^2 > are the centers of these semispheroids. 



For the smooth neck, the following 4th-order expression |]T0| was adopted, namely 

Keck(2;) = ^^eii^li{z - ZifQ{\z\ - \Zi\) , (37) 

where d{x) = for x > and Q{x) = 1 for x < and = —Xjlzl. 
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The frequency Upi in Eq. (^) must be 2;-dependent in order to interpolate smoothly 
between the values uj°^ of the lateral frequencies associated with the left {i = 1) and right 
{i = 2) semispheroids, which are not equal in asymmetric cases. The frequencies 00°^ {i = 
1, 2) characterize the latteral harmonic potentials associated with the two semispheroids 
outside the neck region. In the implementation of such an interpolation, we folUow Ref. 

The angular-momentum dependent term U{lf), where li and I2 are pseudoangular mo- 
menta with respect to the left and right centers zi and Z2, is a direct generalization of the 
corresponding term familiar from the one-center Nilsson potential (e.g., see Ref. [40a]). Its 
function is to lift the usual harmonic-oscillator degeneracies for different angular momenta, 
that is, for a spherical shape the Id — 2s degeneracy is properly lifted into a Id shell that is 
lower than the 2s shell (for the parameters entering into this term, which ensure a proper 
transition from the case of the fissioning cluster to that of the separated two fragments, we 
have followed Ref. [^). 

The cluster shapes associated with the spatial-coordinate-dependent single-particle po- 
tential V{p, z) in the hamiltonian (^) (i.e., the second, thrid, and fourth terms) are deter- 
mined by the assumption that the cluster surface coincides with an equipotential surface of 
value Vo, namely, from the relation V{p, z) = Vq. Subsequently, one solves for p and derives 
the cluster shape p = p{z). For the proper value of Vq, we take the one associated with a 
spherical shape containing the same number of atoms, A^, as the parent cluster, namely, 
Vo = ^meCJo-R^, where fiooo = 4:9r~'^N~^^^ eV, R = TgN^^^, and is the Wigner-Seitz 
radius in atomic units (monovalent metals have been assumed). Volume conservation is 
implemented by requiring that the volume enclosed by the fissioning cluster surface (even 
after separation) remains equal to AnR^/S. 

The cluster shape in this parametrization is specified by four independent parameters. 
We take them to be: the separation d = Z2 — zi of the semispheroids; the asymmetry 
ratio Qas = ^^p2/^pv ^he deformation ratios for the left (1) and right (2) semispheroids 
Qi = uJzi/i^°pi [i = 1,2). 

The single-particle levels of the hamiltonian in Eq. (^) are obtained by numerical 
diagonalization in a basis consisting of the eigenstates of the following auxiliary hamiltonian: 

1 1 

Ho = T+ -meujlp^ + -meujli{z - Zif , (38) 

where uJp is the arithmetic average of and uj°p2- The eigenvalue problem specified by the 
auxiliary hamiltonian (|38D is separable in the cylindrical variables p and z. The general 
solutions in p are those of a two-dimensional oscillator, while in z they can be expressed 



through the parabolic cylinder functions ||72[ . The matching conditions at 2; = for the left 
and right domains yield the 2;-eigenvalues and the associated eigenfunctions [|69[] . 

Having obtained the single-particle spectra, the empirical shell correction (in the spirit 
of Strutinsky's method ||38|), AE'fj^'', is determined from Eq. (|T^) . 

The single-particle average, [i.e., E^p in Eq- (P^)], is calculated |73[ through an 
h expansion of the semiclassical partition function introduced by Wigner and Kirkwood 
(see references in Ref. |7^). For general-shape potentials, this last method amounts ||73| to 



eliminating the semiclassical Fermi energy A from the set of the following two equations 
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(A - vy/^V^V - V{X - \/)-i/2y2^ 



(40) 



where Ne is the total number of delocahzed valence electrons, and V{p, z) is the potential 
in the single-particle hamiltonian of Eq. (|S^). The domain of integration is demarcated by 
the classical turning point r^^, such that VijCy) = A. 

Finally, from the liquid-drop-model contributions, we retain the two most important 
ones, namely the surface contribution and the Coulomb repulsion. To determine the sur- 
face contribution, we calculate numerically the area of the surface of the fissioning cluster 
shape, p = p{z), and multiply it by a surface-tension coefficient specified via an ETF-LDA 
calculation for spherical jellia [^,^] (or even from experimental values). The Coulomb 
repulsion is calculated numerically using the assumption of a classical conductor, namely 
the excess 2 units of positive charge are assumed to be distributed over the surface of the 
fissioning cluster, and in addition each of the fragments carries one unit of charge upon sep- 
aration (for a more elaborate application of the LDM to triaxially deformed ground states 
of neutral and charged metal clusters described via a one-center shape parametrization, see 
our discussion in section II. B. 2 in connection with Eqs. ([T7| - p6D and Ref. [0)- 

As a result, the total energy -Etotai for a specific fission configuration is given by 

-Etotal = -Eldm + '^Ef^^ = Es + Ec + AEg^h'' , 
where Es and Ec are the surface and Coulomb terms, respectively. 



(41) 



III. EXPERIMENTAL TRENDS AND THEORETICAL INTERPRETATION 



In the following, we describe applications of the SE-SCM approach to systematic inves- 
tigations of the effects of shape deformations on the energetics of fragmentation processes 
of metal clusters 



2e 



metal clusters 43 



, and to studies of deformations and barriers in fission of charged 
We mention that, in addition, Strutinsky calculations using phenomeno- 
logical potentials have been reported for the case of neutral sodium clusters assuming axial 
symmetry in Refs. [^^[74| , and for the case of fission in Refs. p8| , ^ |. 



A. Electronic shell effects in monomer and dimer separation energies 



Monomer and dimer separation energies associated with the unimolecular reactions 
— i> K^_^+K, — i> K^_2+K2, and Na^ — >• Na^_;^+Na can be calculated as follows 



N 



E, 



total 



'Z = +l,N -l) + E, 



total 



l)-Etotal(2 = +l,iV) 



(42) 
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FIG. 2. Monomer separation energies, -Dj^^ [see Eq. (|42|)1, from singly cationic clusters 
in the range 5 < < 27. Solid dots: Theoretical results derived from the SE-SCM method. Open 
squares: Experimental measurements from Ref. |75|. Top panel: The spherical model compared 
to experimental data. Middle panel: The spheroidal (axially symmetric) model compared to 
experimental data. Lower panel: The ellipsoidal (triaxial) model compared to experimental data. 



and 

Dt,N = Etot.i{Z = +l,N-2) + Etotai(2 = 0, iV = 2) - Etotai(2 = +1, N) , (43) 

where Z = ±Z {Z being the excess positive or negative charge in absolute units). 

The theoretical results for -D^^ and -D^^r for potassium are displayed in Fig. 2 and Fig. 
3, respectively, and are compared to the experimental measurements [|73|. The theoretical 
and experimental [|76| results for Df^ in the case of sodium are displayed in Fig. 4 (bottom 
panel). An inspection of all three figures leads to the same conclusion as in the case 
of IPs and electron affinities (see Fig. 1 and Ref. [40a]), i.e., that results obtained via 
calculations restricted to spherical shapes compare rather poorly with the experiment, that 
improvement is evident when spheroidal (axially symmetric) deformations are considered, 
and that the agreement between theory and experiment becomes detailed when triaxiality 
(i.e., ellipsoidal shapes) is taken into consideration. The feature of the appearance of strong 
odd-even alternations for iV = 12 — 15 together with a well-defined quartet in the range 
iV = 16 — 19 is present in the experimental monomer separation energies of both potassium 
and sodium clusters, and theoretically it can be accounted for only after the inclusion of 
triaxial deformations. 
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FIG. 3. Dimer separation energies, -D^^v [see Eq. (|43|)], from singly cationic clusters in 
the range 5 < N < 27. Solid dots: Theoretical results derived from the SE-SCM method. Open 
squares: Experimental measurements from Ref. [^. Top panel: The spherical model compared to 
experimental data. Middle panel: The spheroidal model compared to experimental data. Lower 
panel: The ellipsoidal model compared to experimental data. 

We note that in the case of dimer separation energies (Fig. 3) the odd-even alternations 
cancel out. Parents with closed shells or subshells correspond to maxima, while daughters 
with closed shells or subshells are associated with minima (e.g., the triplets = 9 — 11, or 
iV = 15-17). 

We also include for comparison results obtained by KS-LDA calculations [|77| for de- 
formed NaAT clusters restricted to spheroidal (axial) symmetry (Fig. 4, top panel). As 
expected, except for very small clusters (A^ < 9), these results do not exhibit odd-even 
oscillations. In addition, significant discrepancies between the calculated and experimen- 
tal results are evident, particularly pertaining to the amplitude of oscillations at shell and 
subshell closures. 




B. Electronic shell effects in fission energetics 



M^, has attracted considerable attention in 



Fission of doubly charged metal clusters 
the last few years. LDA calculations for fission energetics have usually been restricted 
to spherical jellia for both parent and daughters, |^,19| with the exception of molecular- 
dynamical calculations for sodium [JTOf and potassium [O] clusters with A^ < 12. We 
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FIG. 4. Monomer separation energies, [see Eq. (|4^)], from singly cationic Na^ clusters 

in the range 5 < < 39. Open squares: Experimental measurements from Ref. ||76|. Solid dots 
(Bottom panel): Theoretical results derived from the SE-SCM method in the case of triaxial 
deformations. Solid squares (Top panel): Theoretical results according to the KS-LDA spheroidal 
calculations of Ref . [77|. 



present here systematic calculations for the dissociation energies Aj^^p of the fission pro- 
cesses — i> Kp + K^_p, as a function of the fission channels P. 
We have calculated the dissociation energies 

An,P = EtotA^ = +1, P) + ^total(2 = +1, iV - P) - ^total(2 = +2, A^) , (44) 

for the cases of parent clusters having N = 26, 23, 18, and 15 potassium atoms, and 
compared them with experimental results [|79|. The theoretical calculations compared to 
the experimental results are displayed in Figs. (5 — 8) for = 26, 23, 18, 15, respectively. 
Again, while consideration of spheroidal shapes improves greatly the agreement between 
theory and experiment over the spherical model, fully detailed correspondence is achieved 
only upon allowing for triaxial-shape deformations (notice the improvement in the range 
P = 12- 14 for = 26, and in the range P = 10- 13 for A^ = 23). In the cases A^ = 18 and 
A^ = 15 (Fig. 7 and Fig. 8), the biaxial and triaxial results are essentially identical, since no 
fragment with more than nine electrons is involved. We note that the magic fragments K3" 
and Kg" correspond always to strong minima, and that for A^ = 18 the channel associated 
with the double magic fragments (Kg", Kg") is clearly the favored one over the other magic 
channel with K^, in agreement with the experimental analysis. 

Finally, we carried out calculations of dissociation energies, A^"*^ and A^^^, of the most 
favored fission channels over the whole range up to A^ = 100 atoms for the cases of doubly 
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FIG. 5. Fission dissociation energies, A26,p [see Eq. (^4[)], for the doubly cationic K^q 
cluster as a function of the fission channels P. Solid dots: Theoretical results derived from the 
SE-SCM method. Open squares: Experimental measurements from Ref. [79|. Top panel: The 
spherical model compared to experimental data. Middle panel: The spheroidal model compared 
to experimental data. Lower panel: The ellipsoidal model compared to experimental data. 



charged cationic and anionic sodium clusters, respectively. The triaxial results compared 
to the spherical-jellia calculations according to the LDA-SCM method [^] are displayed 
in Fig. 9 and Fig. 10. In both cases, the main difference from the spherical jellium is 
a strong suppression of the local minima, indicating that the critical size for exothermic 
fission is significantly smaller than = 100 (about = 30), as indeed has been observed 
experimentally for hot cationic alkali-metal clusters (the spherical-jellium results clearly 
are not compatible with the emergence of such experimental critical sizes in the size range 
< 100). 



C. Electronic shell effects in fission barriers and fission dynamics of metal clusters 

In this section, we focus our discussion on recent trends in studies of binary fission 
processes in doubly charged metal clusters. 
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FIG. 6. Fission dissociation energies, A23,p [see Eq. (^4|)], for the doubly cationic K^^ 
cluster as a function of the fission channels P. Solid dots: Theoretical results derived from the 
SE-SCM method. Open squares: Experimental measurements from Ref. [79|. Top panel: The 
spherical model compared to experimental data. Middle panel: The spheroidal model compared 
to experimental data. Lower panel: The ellipsoidal model compared to experimental data. 



1. Molecular- dynamics studies of fission 



Before discussing applications of the SE-SCM (and variants thereof) to the description 
of cluster fission, we note that for atomic and molecular clusters microscopic descriptions of 
energetics and dynamics of fission processes, based on modern electronic structure calcula- 
tions in conjunction with molecular dynamics simulations (where the classical trajectories of 
the ions, moving on the concurrently calculated Born-Oppenheimer (BO) electronic poten- 
tial energy surface, are obtained via integration of the Newtonian equations of motion), are 
possible and have been performed [0,0 using the BO-local-spin-density(-LSD)-functional- 
MD method |]80|. Such calculations, using norm- conserving non-local pseudopotentials 
and self-consistent solutions of the KS-LSD equations [|l3,0, applied to small sodium 
TI| and potassium |rT| clusters, revealed several important trends (Figs. 11 — 13): (i) 



The energetically favorable fission channel for such doubly-charged clusters is the asym- 
M]^_3+ M^, containing a "magic" daughter (M = Na, K), i.e., 



metric one 



N > Nl 



where N^'^ 



p) + EiMj>) - E{W^) is smallest for P = 3; (ii) Fission of clusters with 
7, involves barriers, whose magnitudes reflect the closed-shell 
stability of the parent cluster (i.e., Ei, for = 10 is particularly high), exhibiting a double- 
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FIG. 7. Fission dissociation energies, Aig^p [see Eq. (^4[)], for the doubly cationic K^^ 
cluster as a function of the fission channels P. Solid dots: Theoretical results derived from the 



SE-SCM method. Open squares: Experimental measurements from Ref. [79|. Top panel: The 
spherical model compared to experimental data. Middle panel: The spheroidal model compared 
to experimental data. Lower panel: The ellipsoidal model compared to experimental data. 



humped barrier shape [see, Figs. 11 and 13(a)]; (ill) The eventual fission products may be 
distinguishable (i.e., preformed) already at a rather early stage of the fission process (on 



the top of the exit barrier for Na^^Q , see Fig. 12, or prior to the exit barrier for see 
Fig. 13), and the electronic binding between the two fragments is long-range in nature; (iv) 
The kinetic energy release Sr in the favorable channel obtained via dynamic simulations 
was found to be given by £r ^ Eb + \ An^3\, and the results are in correspondence with 
experimental measurements |lll] for (5 < < 12). Furthermore, in agreement with 
experimental findings, the emerging fragments are vibrationally excited, with the heating 
of the internal nuclear degrees of freedom of the fission products in the exit channel origi- 
nating from dynamical conversion of potential into internal kinetic energy [see. Kg",!: in Fig. 
13(b)]. 



2. SE-SCM interpretation of fissioning processes 

The method we adopt in this section for further studying metal-cluster fission is the 
SE-SCM described in section II. B. 6 (see also Ref. [Q). 
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FIG. 8. Fission dissociation energies, Ais^p [see Eq. (^4|)], for the doubly cationic Kf^ 
cluster as a function of the fission channels P. Solid dots: Theoretical results derived from the 



SE-SCM method. Open squares: Experimental measurements from Ref. [79|. Top panel: The 
spherical model compared to experimental data. Middle panel: The spheroidal model compared 
to experimental data. Lower panel: The ellipsoidal model compared to experimental data. 



As discussed above (see section II. B), in the SE-SCM method we need to introduce 
appropriate empirical potentials. As will become apparent from our results, one-center 
potentials (like the one-center modified, anisotropic harmonic oscillator) are not adequate 
for describing shell effects in the fission of small metal clusters; rather, a two-center potential 
is required. Indeed, the empirical potentials should be able to simulate the fragmentation 
of the initial parent cluster towards a variety of asymptotic daughter-cluster shapes, e.g., 
two spheres in the case of double magic fragments, a sphere and a spheroid in the case of a 
single magic fragment, or two spheroids in a more general case. In the case of metal clusters, 
asymmetric channels are most common, and thus a meaningful and flexible description of 
the asymmetry is of primary concern. We found |^3[ that such a required degree of flexibility 
can be provided via the shape parametrization of the asymmetric two-center-oscillator shell 
model (ATCOSM) introduced earlier in nuclear fission [39|| (see section II. B. 6). 

In addition to the present shape parametrization other two-center shape 

parametrizations [mainly in connection with KS-LDA jellium calculations] have been used 
||8T|^3| in studies of metal cluster fission. They can be grouped into two categories, namely, 
the two-intersected-spheres jellium |ET|,p3|, and the variable-necking-in parametrizations 



82,83]. In the latter group, Ref. [82] accounts for various necking-in situations by using 
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FIG. 9. Solid dots: LDA-SCM results for the dissociation energies A^"^ for the most 
favorable fission channel for doubly charged cationic parents Na^ when the spherical jellium is 
used. The influence of triaxial deformation effects (calculated with the SE-SCM approach) is 
shown by the thick dashed line. 
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FIG. 10. Solid dots: LDA-SCM results for the dissociation energies for the most 

favorable fission channel for doubly charged anionic parents Na^ when the spherical jellium is 
used. The influence of triaxial deformation effects (calculated with the SE-SCM approach) is 
shown by the thick dashed line. 



the "funny-hills" parametrization 



183| describes the necking-in by con- 



while Ref. 

necting two spheres smoothly through a quadratic surface. The limitation of these other 
parametrizations is that they are not flexible enough to account for the majority of the 
effects generated by the shell structure of the parent and daughters, which in general do 
not have spherical, but deformed (independently from each other), shapes. An example is 
offered by the case of the parent Na^g", which has a metastable oblate ground state, and 
thus cannot be described by any one of the above parametrizations. We wish to emphasize 
again that one of the conclusions of the present work is that the shell structures of the 
(independently deformed) parent and daughters are the dominant factors specifying the 



fission barriers, and thus parametrizations |]81|-p^ with restricted final fragment (or par- 
ent) shapes are deficient in accounting for some of the most important features governing 
metal-cluster fission. 

As a demonstration of our method, we present results for two different parents, namely 
Na?j[ and Na?^. 
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FIG. 11. Molecular dynamics results for the potential energy vs distance (in atomic units) 
between the centers of mass for the fragmentation of Na^(j" into Na^ and Na^ (sohd) and Na^ 
and Na"*" (dashed), obtained via constrained minimization of the LSD ground-state energy of the 
system ||To[| . 



Fig. 14 presents results for the channel Na^J" — > Nay + Na;^ for three different cases, 
namely, when the larger fragment Na^" is oblate (left column), spherical (middle column), 
and prolate (right column). From our one-center analysis, we find as expected that Na^" 
(with six electrons) has an oblate global minimum and a higher in energy prolate local 
minimum. In the two-center analysis, we have calculated the fission pathways so that the 
emerging fragments correspond to possible deformed one-center minima. It is apparent that 
the most favored channel (i.e., having the lowest barrier, see the solid line in the bottom 
panels) will yield an oblate Na^ (left column in Fig. 14), in agreement with the expectations 
from the one-center energetics analysis. 

The middle panels exhibit the decomposition of the total barrier into the three compo- 
nents of surface. Coulomb, and shell-correction terms [see Eq. (|4lD ], which are denoted by 
an upper dashed curve, a lower dashed curve, and a solid line, respectively. The total LDM 
contribution (surface plus Coulomb) is also exhibited at the bottom panels (dashed lines). 

It can be seen that the LDM barrier is either absent or very small, and that the total 
barrier is due almost exclusively to electronic shell effects. The total barrier has a double- 
humped structure, with the outer hump corresponding to the LDM saddle point, which 
also happens to be the scission point (indicated by an empty vertical arrow). The inner 
hump coincides with the peak of the shell-effect term, and is associated with the rearrange- 
ment of single-particle levels from the initial spherical parent to a molecular configuration 
resembling a Nay attached to a Naj^. Such molecular configurations (discovered earlier 
in first-principles MD simulations |TO| , p]l|] of fission of charged metal clusters, as well as 
in studies of fusion of neutral clusters [^) are a natural precursor towards full fragment 
separation and complete fission, and naturally they give rise to the notion of preformation 



of the emerging fragments [|TO| , pT 



Fig. 15(a) displays the ATCOSM results for the symmetric channel Na^^ — >■ 2Nag, when, 
for illustrative purposes, the parent is assumed to be spherical at d = (this channel is 
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FIG. 12. Pragmentation dynamics of Na^^ from first-principles Born-Oppenheimer Lo- 
cal-spin-density functional Molecular-dynamics simulations |10|. (a)-(c) Center-of-mass distance 



between the eventual fission products (-Rc.m.-c.m.)) total potential energy {Ep), and the electronic 
contribution Eg to Ep, vs time, (d)-(f) Contours of the total electronic charge distribution at 
selected times calculated in the plane containing the two centers of mass. The axis is parallel 
to Rc.m.-c.m.- (g)-(i) Cluster configurations for the times given in (d)-(f). Dark and light balls 
represent ions in the large and small fragments, respectively. Energy, distance, and time in units 
of eV, bohr (oq), and ps, respectively. 



favored compared to that of the trimer , both from energetics and barrier considerations; 
for small clusters, this is the only case where a channel other than that of the trimer is 
the most favored one). The top panel of Fig. 15(a) describes the evolution of the single- 
particle spectra. The spherical ordering Is, Ip, Id, 2s, etc., for the parent at d = is 
clearly discernible. With increasing separation distance, the levels exhibit several crossings, 
and, after the scission point, they naturally regroup to a new ordering associated with the 
spherical Nag products (at the end of the fission process, the levels are doubly degenerate 
compared to the initial configuration, since there are two Nag" fragments). It is seen that 
the ATCOSM leads to an oscillator energy (i.e., the gap between two populated major shells 
exhibited at the right end of the figure) of 1.47 eV for each Nag fragment in agreement 
with the value expected from the one-center model [the Is state of Nag" lies at 2.21 eV; in 
the case of the initial spherical Na^g" {d = 0), the oscillator energy corresponding to the 
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FIG. 13. (a) Potential energy of K.I2 fissioning in tlie favorable channel (Kjj'+K^) versus 
the inter-fragment distance -Rc.m.-c.m. obtained via constrained minimization. The origin of the 
Ep scale is set at the optimal pre-barrier configuration (A). For large i?c.m.-c.m., Ep=— 0.9 eV, 
i.e., A3. Included also are cluster configurations of K^g^ corresponding to: a compact isomer (A') 
(the energy of the optimal compact isomer found is denoted by an arrow); the optimal bound 
configuration (A) ; the structure on top of the exit-channel barrier for which contours of the total 
electronic charge density, p, are shown yil[ . 

(b) Time evolution of i?c.m.-c.m.5 the internal vibrational kinetic energies of the fragments (K^^^ 
and K|j°+*) and the sum of the fragments translational kinetic energies (i^^cm) obtained via a 
BO-LSD-MD simulation starting from ionization (t = 0) of a cluster at 500 K. A line is 
drawn in KJjJ (for t > 3 ps) to guide the eye, illustrating heating of the internal vibrational 
degrees of freedom of the departing fragment. 
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FIG. 14. ATCOSM results for the asymmetrie channel Na^J Naif +Na^. The final con- 
figuration of Nag' is spherical. For the heavier fragment Na^", we present results associated with 
three different final shape configurations, namely, oblate [(o,s); left], spherical [(s,s); middle], and 
prolate [(p,s); right]. The ratio of shorter over longer axis is 0.555 for the oblate case and 0.75 for 
the prolate case. 

Bottom panel: LDM energy (surface plus Coulomb, dashed curve) and total potential energy 

(LDM plus shell corrections, solid curve) as a function of fragment separation d. The empty 
vertical arrow marks the scission point. The zero of energy is taken at d = 0. A number ( — 1.58 
eV or —0.98 eV), or a horizontal solid arrow, denotes the corresponding dissociation energy. 
Middle panel: Shell-correction contribution (solid curve), surface contribution (upper dashed 
curve), and Coulomb contribution (lower dashed curve) to the total energy, as a function of 
fragment separation d. 

Top panel: Single-particle spectra as a function of fragment separation d. The occupied (fully or 
partially) levels are denoted with solid lines. The unoccupied levels are denoted with dashed lines. 
On top of the figure, four snapshots of the evolving cluster shapes are displayed. The solid ver- 
tical arrows mark the corresponding fragment separations. Observe that the doorway molecular 
configurations correspond to the second snapshot from the left. Notice the change in energy scale 
for the middle and bottom panels, as one passes from (o,s) to (s,s) and (p,s) final configurations. 
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FIG. 15. ATCOSM results for the symmetric channel Na^g" — ^ 2Nag", when the initial 
parent shape is assumed (a) spherical, and (b) oblate (with a shorter over longer axis ratio equal 
to 0.699). Panel distribution and other notations and conventions are the same as in Fig. 14. 
The top dotted line in the bottom panel of (a) represents the total energy without the Coulomb 
contribution. Observe that the doorway molecular configurations correspond to the third snapshot 
from the left. Notice that the zero of all energies is taken at d = 0. 



gap between major shells is 1.17 eV, and the corresponding Is state lies at 1.75 eV]. 

From the middle panel of Fig. 15(a), we observe that the shell-correction (solid line) 
contributes a net gain in energy of about 1.6 eV upon dissociation into two Nag" fragments. 
This gain is larger than the increase in energy (i.e., positive energy change) due to the surface 
term, which saturates at a value of about 1 eV after the scission point at (i ~ 23 a.u. The 
total energy is displayed in the bottom panel of Fig. 15(a) (solid line) along with the LDM 
barrier (dashed line). Even though distorted (when compared to the cases of Fig. 14), the 
total barrier still exhibits a two-peak structure, the inner peak arising from the hump in the 
shell correction, and the outer peak arising from the point of saturation of the surface term 
(this last point coincides again with the scission point, as well as with the saddle of the LDM 
barrier). An inner local minimum is located at d ^ 8 a.u., and corresponds to a compact 
prolate shape of the parent [see second drawing from the left at the top of Fig. 15(a)], while a 
second deeper minimum appears at d ~ 18 a.u., corresponding to a superdeformed shape of 
a molecular configuration of two Nag" clusters tied up together [preformation of fragments, 
see third drawing from the left at the top of Fig. 15(a)]. The inner barrier separating 
the compact prolate configuration from the superdeformed molecular configuration arises 
from the rearrangement of the single-particle levels during the transition from the initially 
assumed spherical Na^g configuration to that of the supermolecule Nag -l-Nag . We note 
that the barrier separating the molecular configuration from complete fission is very weak 
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being less than 0.1 eV. 

The top dotted hne at the bottom panel displays the total energy in the case when the 
Coulomb contribution is neglected. This curve mimics the total energy for the fusion of 
two neutral Nag clusters, namely the total energy for the reaction 2Na8 Naig. Overall, 
we find good agreement with a KS-LDA calculation for this fusion process (see Fig. 1 of 
Ref. [^). We further note that the superdeformed minimum for the neutral Naie cluster 
is deeper than that in the case of the doubly charged Na^^ cluster. Naturally, this is due 
to the absence of the Coulomb term. 

The natural way for producing experimentally the metastable Na^g" cluster is by ioniza- 
tion of the stable singly-charged Na^g cluster. Since this latter cluster contains seventeen 
electrons and has a deformed oblate ground state [40a], it is not likely that the initial con- 
figuration of Na^^ will be spherical or prolate as was assumed for illustration purposes in 
Fig. 15(a). Most likely, the initial configuration for Na^^ will be that of the oblate Na^'g. 
To study the effect that such an oblate initial configuration has on the fission barrier, we 
display in Fig. 15(b) ATCOSM results for the pathway for the symmetric fission channel, 
starting from an oblate shape of Na^^, proceeding to a compact prolate shape, and then 
to full separation between the fragments via a superdeformed molecular configuration. We 
observe that additional potential humps (in the range 2 a.u. < d < 6 a.u.), associated with 
the shape transition from the oblate to the compact prolate shape, do develop. Concerning 
the total energies, the additional innermost humps result in the emergence of a significant 
fission barrier of about 0.52 eV for the favored symmetric channel [see d ^ 5 a.u. in Fig. 
15(b)]. 

From the above analysis, we conclude that considerations of the energy pathways lead- 
ing from the parent to preformation configurations (i.e., the inner-barrier hump, or humps) 
together with the subsequent separation processes are most important for proper elucida- 
tion of the mechanisms of metal-cluster fission processes. This corroborates earlier results 
obtained via first-principles MD simulations [ pD| , pT| pertaining to the energetics and dynam- 
ical evolution of fission processes, and emphasizes that focusing exclusively |pT| , |83| on the 
separation process between the preformed state and the ultimate fission products provides 
a rather incomplete description of fission phenomena in metal clusters. It is anticipated 
that, with the use of emerging fast spectroscopies 1^^, experimental probing of the detailed 
dynamics of such fission processes could be achieved. 



IV. INFLUENCE OF ELECTRONIC ENTROPY ON SHELL EFFECTS 

In the previous sections, we showed that consideration of triaxial (ellipsoidal) shapes in 
the framework of the SCM leads to overall substantial systematic improvement in the agree- 
ment between theory and experimental observations pertaining to the major and the fine 
structure of the size-evolutionary patterns associated with the energetics of fragmentation 
processes ( monomer /dimer dissociation energies and fission energetics) and ionization. 

The theoretical methods and discussion of deformation effects in the previous sections 
were restricted to zero temperature. However, the experiments are necessarily made with 
clusters at finite temperatures, a fact that strongly motivates the development of finite- 
temperature theoretical approaches. 
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FIG. 16. Monomer separation energies of clusters at two temperatures, T = 10 K, and 
300 K. Solid dots: Theoretical FT-SE-SCM results. Open squares: experimental measurements 
1 )7 5| . To facilitate comparison, the SE-SCM results at the higher temperature have been shifted 
by 0.07 eV, so that the theoretical curves at both temperatures refer to the same point at = 10. 



Due to the difficulty of the subject, to date only a few finite-temperature theoretical 
studies of metal clusters have been performed. In this section, we discuss briefiy some of 
the conclusions of a recent SCM study regarding the importance of thermal effects. 



The theoretical details pertaining to this finite-temperature (FT) — SE— SCM will not be 
elaborated here, but they can be found in the aforementioned reference. 

The main conclusion of Ref. [0] was that, in conjunction with deformation effects, 
electronic-entropy effects in the size-evolutionary patterns of relatively small (as small as 
20 atoms) simple-metal clusters become prominent already at moderate temperatures. At 
smaller sizes, electronic-entropy effects are less prominent, but they can still be discernible. 
As an example, we present in Fig. 16 the monomer separation energies of clusters 
for two temperatures (T = 10 K and T = 300 K), along with the available experimental 
measurements (open squares) in the size range = 4 — 23. First notice that the 
T = 10 K results are practically indistinguishable from the T = K results presented in 
Fig. 2. Compared to the T = 10 K results, the theoretical results at T = 300 K are in 
better agreement with the experimental ones due to an attenuation of the amplitude of 
the alternations (e.g., notice the favorable reduction in the size of the drops at iV = 9, 15, 
and 21). This amplitude attenuation, however, is moderate, and it is remarkable that the 
T = 300 K SCM results in this size range preserve in detail the same relative pattern as the 
T = K ones (in particular, the well-defined odd-even oscillations in the range = 4 — 15 
and the ascending quartet at iV = 16 — 19 followed by a dip at = 20). 

As a further example, the theoretical IPs of Kat clusters in the size range 3 < A^ < 102 for 
three temperatures, T = 10 K, 300 K, and 500 K, are displayed in Fig. 17, and are compared 
with the experimental measurements [^] (open squares; the experimental uncertainties are 
0.06 eV for A^ < 30 and 0.03 eV for A^ > 30). As was the case with our earlier T = 



K results P0[, the T = 10 K theoretical results exhibit the following two characteristics: 



(i) Above A^ = 21, a pronounced fine structure between major-shell closures which is not 
present in the experimental measurements; (ii) Steps at the major-shell closures which are 
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FIG. 17. IPs of Ktv clusters at three temperatures, T = 10 K, 300 K, and 500 K. Solid dots: 



Theoretical FT-SE-SCM results. Open squares: experimental measurements [36|. 



much larger than the experimental ones, i.e., three-to-five times for = 40, 58, and 92, 
and two-to-three times for = 8 and 20 [this needs to be contrasted to the experimental 
IPs for cold Na^v clusters, which are in overall good agreement with out T = K SE-SCM 
results regarding both characteristics (see Fig. 1)]. 

The agreement between theory and experiment is significantly improved at T = 300 
K. Indeed, in comparison with the lower-temperature calculations, the T = 300 K results 
exhibit the following remarkable changes: (i) Above A^ = 21, the previously sharp fine- 
structure features are smeared out, and as a result, the theoretical curve follows closely 
the mild modulations of the experimental profile. In the size range A^ = 21 — 34, three 
rounded, hump-like formations (ending to the right at the subshell closures at A^ = 26, 
30, and 34) survive in very good agreement with the experiment (the sizes of the drops at 
A^ = 26, 30 and 34 are comparable to the experimental ones [^); (ii) The sizes of the IP 
drops at A^ = 20, 40, 58, and 92 are reduced drastically and are now comparable to the 
experimental ones. In the size range A^ < 20, the modifications are not as dramatic. Indeed, 
one can clearly see that the pattern of odd-even alternations remains well defined, but with 
a moderate attenuation in amplitude, again in excellent agreement with the experimental 
observation. 

For T = 500 K, the smearing out of the shell structure associated with the calculated 
results progresses even further, obliterating the agreement between theory and experiment. 
Specifically, the steps at the subshell closures at A^ = 26 and 30, as well as at the major- 
shell closures at A^ = 40, 58, and 92 are rounded and smeared out over several clusters (an 
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analogous behavior has been observed in the logarithmic abundance spectra of hot, singly 
cationic, copper, silver, and gold clusters |^). At the same time, however, the odd-even 
alternation remains well defined for < 8. We further notice that, while some residue 
of fine structure survives in the range = 9 — 15, the odd-even alternations there are 
essentially absent (certain experimental measurements ||90| of the IPs of hot NaAr clusters 
appear to conform to this trend). 

The infiuence of the electronic entropy on the height of fission barriers has not been 
studied as yet, but it will undoubtedly be the subject of future research in metal-cluster 
physics. In any case, based on the results of this section, it is natural to conjecture that 
electronic-entropy effects will tend to quench the barrier heights, especially in the case of 
larger multiply charged clusters. 



V. SUMMARY 

In this chapter, we have elucidated certain issues pertaining to evaporation and fission 
processes of metallic clusters, focusing on electronic shell effects and their importance in 
determining the energetics, structure, pathways, and dynamical mechanisms of dissociation 
and fragmentation in these systems, and have outlined and demonstrated various theoretical 
approaches currently used in investigations of cluster fragmentation phenomena, ranging 
from microscopic first-principles electronic structure calculations coupled with molecular 
dynamics simulations to adaptation of more phenomenological in nature models originated 
in studies of atomic nuclei. In this respect, a recurrent theme in this exposition has been 
the crucial importance of deformation and electronic-entropy (temperature) effects, as well 
as their treatment with the help of shell correction methods. 

By drawing analogies, as well as differences, between certain aspects of nuclear fission 
and nuclear radioactivity phenomena and atomic (metallic) cluster fission processes, we 
have attempted to provide a unifying conceptual framework for discussion of the physi- 
cal principles underlying modes of cluster fission (i.e., importance of deformations, shell 
effects originating from fragments and parent, asymmetric and symmetric fission, single 
and double-humped barriers, fissioning cluster shapes, and dynamical aspects, such as the 
time-scale of fission processes, kinetic energy release, and dynamical energy redistribution 
among the fission products). 

We conclude by commenting on some experimental and theoretical issues in cluster 
fission which remain as future challenges (limiting ourselves to metallic clusters). These 
include: fission dynamics of multiply charged large metal clusters []6|, p3| , |2^ ; systematic in- 
vestigations of temperature effects on modes of cluster fission, and ternary, and higher 
multi-fragmentation processes; time-resolved spectroscopy of fission processes and of fission 
isomers; spin effects in fission; tunneling processes and corresponding life-times in sub- 
barrier fission modes of clusters of light elements, e.g., lithium; and fission processes of 
non-simple metal clusters. 
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